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Abstract.  We  study  the  control  of  the  rotation  of  a laser-cooled  ion  crystal  in  a Penning  trap  by 
a rotating  electric  field  perturbation.  We  show  that  application  of  a small  torque  produces  sudden 
angular  jumps  or  ‘slips’  of  the  crystal  orientation  spaced  by  intervals  when  the  crystal  is  phase- 
locked  or  ‘stuck’  relative  to  the  rotating  perturbation.  The  distribution  of  angular  slips  is  described 
by  a power  law,  where  the  power-law  exponent  depends  on  the  applied  torque.  We  believe  this 
system  is  driven  by  a constant  force  and  that  small  perturbations  or  thermal  effects  trigger  the  slips. 


Non-neutral  plasmas  confined  in  Penning-Malmberg  traps  are  used  in  a variety  of 
experiments  including  plasma  physics,  Coulomb  crystal  studies,  precision  spectroscopy, 
antimatter  research,  and  storage  of  highly  charged  ions  [1],  Recently  there  has  been 
a great  deal  of  interest  in  using  a rotating  electric  field  perturbation  to  control  the 
global  E x B rotation  of  these  plasmas  [2,  3,  4].  For  crystallized  ion  plasmas,  phased- 
locked  control  of  the  plasma  rotation  has  been  demonstrated  [4],  which  has  important 
implications  for  atomic  clocks  [5]  and  for  quantum  computation  with  trapped  ions  [6], 

In  this  manuscript  we  study  the  limits  to  phase-locked  control  due  to  the  application 
of  a small  torque  produced  by  the  radiation  pressure  of  a weak  laser  beam.  We  observe 
sudden  angular  jumps  or  ‘slips’  of  the  crystal  orientation  spaced  by  intervals  when 
the  crystal  orientation  is  phase-locked  or  ‘stuck’  relative  to  the  rotating  perturbation. 
Stick-slip  behavior  similar  to  that  observed  here  is  found  in  many  different  and  diverse 
systems:  for  example,  in  studies  of  friction  between  two  surfaces  [7,  8],  in  experiments 
on  avalanches  and  slips  in  granular  systems  [9,  10,  11],  and  as  the  underlying  process 
in  spring-block  models  of  earthquakes  [12,  13].  Many  of  these  systems,  including  the 
study  presented  here,  exhibit  a power-law  distribution  of  the  slip  amplitudes,  indicative 
of  an  underlying  critical  point  [14,  15]. 

Our  work  uses  the  Penning-Malmberg  trap  at  NIST  to  store  ~ 20, 000  9Be+  ions.  The 
ions  interact  by  unscreened  Coulomb  repulsion  and  are  Doppler  laser-cooled  [16]  to 
millikelvin  temperatures,  where  they  become  strongly  coupled.  The  degree  of  coupling 
is  quantified  by  the  Coulomb  coupling  parameter, 
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FIGURE  1.  (a)  Schematic  of  the  cylindrical  Penning  trap  and  the  top-view  imaging  system.  The  side- 
view  imaging  system  is  not  shown,  (b)  Strobed  top-view  image  of  a 5 axial  plane  9Be+  ion  crystal  with  a 
bee  structure,  similar  to  those  used  in  this  study,  (c)  Side-view  image  (unstrobed)  of  the  same  ion  crystal. 
The  diameter  (2/  «,.)  of  the  <)Be+  ions  is  495  fim.  Ions  of  greater  mass  are  located  at  /•  > /•&.  but  do  not 
fluoresce  in  the  laser  beam.  The  rotation  axes  arc  indicated. 


which  is  the  ratio  of  the  Coulomb  potential  energy  of  neighboring  ions  to  the  kinetic 
energy  per  ion.  Here,  £q  is  the  permittivity  of  the  vacuum,  e is  the  charge  of  an  ion,  kg  is 
Boltzmann’s  constant,  T is  the  temperature,  and  am  is  the  Wigner-Seitz  radius,  defined 
by  47t(a„.j)3/3  = 1 //to,  where  no  is  the  ion  number  density.  With  coupling  parameters 
greater  than  T ~ 1 the  ions  are  strongly  coupled  [17],  and  when  T>172  they  form  a 
Coulomb  crystal  (a  classical  Wigner  crystal)  [18,  19],  Structurally  similar  Coulomb 
crystals  are  believed  to  exist  in  dense  astrophysical  matter,  such  as  the  interior  of  white 
dwarfs  and  the  outer  crust  of  neutron  stars  [20].  Interestingly,  observations  of  power-law 
statistics  of  soft  y-ray  (SGR)  events  have  recently  been  interpreted  as  evidence  that  the 
Coulomb  crystal  comprising  the  outer  crust  of  a magnetized  neutron  star  can  undergo 
very  large-scale  slips  (‘starquakes’)  [21,  22],  The  measured  power  law  exponents  of  the 
neutron  starquakes  lie  within  the  range  of  exponents  we  measure  here  (see  Fig.  3). 

Figure  1(a)  shows  the  experimental  setup  [4,  19],  We  create  9Be+  ions  by  ionizing 
9Be  atoms  in  a separate  trap  (not  shown)  and  then  transferring  the  ions  to  the  main  trap 
for  experimentation  [23],  The  9Be+  ions  were  confined  radially  by  a uniform  magnetic 
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field  B=4.465  T (cyclotron  frequency  Qc/2n  = 7.608  MHz)  in  the  z direction  and 
axially  by  a potential  difference  of  Vq  = -500  V applied  between  the  center  and  end 
electrodes  of  the  trap.  Near  the  trap  center  the  trap  potential  is  quadratic  and  given  by 
mco?(z2  — r2 /2)/{2e),  where  the  axial  frequency  to;/27t  = 565  kHz  for  9Be+.  Here  r and 
z denote  the  cylindrical  radius  and  axial  coordinates.  Due  to  the  axial  magnetic  field  and 
the  radial  components  of  the  ion  space  charge  and  trap  electric  fields,  the  ion  crystal 
rotates  at  a frequency  to,  about  the  trap  symmetry  axis  z. 

In  our  work  the  plasma  has  had  sufficient  time  to  evolve  to  a thermal  equilibrium  state 
where  its  rotation  and  density  is  uniform  across  the  cloud  [24],  For  20,000  trapped  ions 
with  ©,./ 2k  = 22.8  kHz  (typical  of  this  work),  the  ion  plasma  has  a density  of  7.1  x 107 
cm-3  and  a lenticular  shape  with  a radial  diameter  2ro  ~ 1 mm  and  an  aspect  ratio 
a = zoAo  ~ 0.05.  Here,  2zo  denotes  the  axial  extent  of  the  plasma  at  r = 0.  In  addition 
to  9Be+  ions,  ions  of  greater  mass  (‘heavy  ions’)  such  as  BeH+  and  BeOH+  are  created 
by  reactions  with  9Be+  ions  and  background  neutral  molecules.  For  the  work  discussed 
here,  typically  20  % to  50  % of  the  plasma  consisted  of  heavy  impurity  ions.  These 
ions  are  sympathetically  cooled  to  temperatures  similar  to  the  9Be+  ions  and,  due  to  the 
rotation,  centrifugally  separate  to  larger  radii,  where  they  crystallize  [25,  26]. 

To  control  to,.,  we  applied  an  electric  field  perturbation  rotating  about  the  z axis  at 
frequency  (£>rp  [4].  The  rotating  perturbation  applies  a torque  on  the  radial  boundary  of 
the  plasma  (on  the  non-fluorescing,  heavy  ions)  by  creating  a small-amplitude  traveling 
wave.  The  torque  due  to  this  wave  is  then  transferred  to  the  plasma  interior  through 
the  strong  inter-particle  Coulomb  forces,  which  act  to  bring  the  plasma  to  the  same 
rotation  frequency  as  (s)rp  [24].  We  observe  similar  stick-slip  motion  with  both  dipole 
and  quadrupole  rotating  fields.  However,  most  measurements,  including  those  we  report 
here,  were  taken  with  a dipole  rotating  field. 

The  radial  binding  force  of  the  trap  is  due  to  the  Lorentz  force  produced  by  the 
plasma’s  rotation  through  the  magnetic  field.  Therefore  changing  (0,  changes  the  radial 
binding  force  of  the  trap  and  provides  a sensitive  way  to  adjust  the  overall  shape  and 
structural  phase  of  the  plasma.  In  this  work,  the  frequency  of  the  rotating  perturbation 
(Orp  was  kept  near  271  x 22.8  kHz  and  slightly  adjusted  to  produce  a disk-shaped  plasma 
consisting  of  5 axial  planes  with  a bcc-like  crystal  structure  in  the  plasma  center  [19], 
Because  ©,■  <C  Qc,  the  ion  motion  in  a direction  perpendicular  to  the  magnetic  field  is 
determined  principally  by  E x B guiding  center  dynamics  [27], 

The  main  cooling-laser  beam  (?c  = 313  nm)  was  directed  along  the  z axis  and  tuned  10 
to  20  MHz  lower  in  frequency  than  a hyperfine-Zeeman  component  of  the  2s  2S]  — ► 2 p 

1Psj2  resonance  with  a natural  linewidth  of  19  MHz.  This  beam’s  power  was  ~50  /rW 
and  the  beam  was  focused  to  a ~0.5  mm  waist  at  the  ion  crystal.  A second  cooling 
beam  (_L  beam  in  Fig.  1(a)),  derived  from  the  same  laser,  was  directed  perpendicularly 
to  z and  had  a ~70  /;m  waist  and  ~1  ji W power.  Both  the  perpendicular  and  parallel 
cooling  lasers  were  required  to  form  a well  defined  crystal  in  the  disk-shaped  plasmas 
discussed  here.  Although  we  did  not  measure  the  ion  temperature,  the  theoretical  cooling 
limit  is  0.5  mK,  and  an  experimental  upper  bound  of  T < 10  mK  has  been  measured 
in  previous  experiments  [23];  for  a density  of  no  = 7.1  x 107  cm-3,  these  limits  give  a 
range  of  1 12  < F < 2230.  An  indirect  estimate  of  the  ions’  coupling  parameter,  based  on 
observations  of  when  the  bcc(l  10)  structural  phase  was  exclusively  favored,  is  T ~ 600 
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[28]. 

The  _L  beam  is  normally  directed  through  the  radial  center  (/-  = 0)  of  the  crystal  in 
order  to  minimize  its  applied  torque  while  providing  a low  Doppler-cooling  temperature 
[16].  In  this  experiment,  we  offset  the  J_-beam  position  slightly  (5  to  30  /im)  from  the 
plasma  center  to  produce  a torque  on  the  9Be+  ions  in  the  same  direction  as  the  plasma 
rotation.  The  resultant  torque  from  the  ± beam  was  larger  than  any  other  ambient  torque 
due  to,  for  example,  resonances  between  modes  and  asymmetries  in  the  trap  construction 
[29]  or  background  gas  drag  [30]. 

A series  of  lenses  formed  side-  and  top-view  images  of  the  ion  fluorescence,  with 
viewing  directions  perpendicular  and  parallel  to  the  magnetic  field  respectively,  on  either 
a gateable  charge-coupled  device  (CCD)  camera,  or  on  an  imaging  photomultiplier 
tube.  The  resolution  of  the  optical  systems  was  ~4  //m,  while  typical  interparticle 
spacings  were  ~15  /mi.  By  detecting  the  ion  fluorescence  synchronously  with  the 
rotating  perturbation  drive,  images  of  the  individual  ions  that  make  up  the  Coulomb 
crystals  were  obtained.  Figure  1(b)  shows  a strobed,  top-view  CCD  camera  image  of  a 
5-axial  plane  crystal  in  the  bcc  structural  phase,  accumulated  over  40  s.  The  ion  positions 
are  well  localized  in  the  plasma  center;  however,  at  larger  radii  they  are  blurred. 

To  investigate  the  blurring  we  used  an  imaging  photomultiplier  tube  in  the  top-view 
position  to  record  the  positions  and  detection  times  of  the  fluorescence  photons.  Data 
sets  consisted  of  125  ms  intervals  of  data  recorded  each  second  over  long  periods  of 
time  (up  to  5000  s).  Images  similar  to  those  in  Fig.  1(b)  were  created  for  each  125 
ms  interval  by  constructing  2D  histograms  of  the  ion  fluorescence  in  the  frame  of  the 
rotating  perturbation.  The  orientation  0m  of  the  central  crystallized  region  in  the  rotating 
frame  was  determined  (modulo  n due  to  the  bcc  crystal  bilateral  symmetry)  with  an 
uncertainty  of  ~ 0.0027t  radians  [31]. 

In  Fig.  2 we  plot  0cr>(t)  for  two  data  sets  that  differ  mainly  in  the  amount  of  -L-beam 
torque.  Over  long  time  scales  the  I-beam  torque  produces  a slightly  faster  rotation  (a 
rotational  ‘creep’)  of  the  9Be+  crystal  relative  to  the  rotating  perturbation.  For  example, 
in  data  set  2 Aco  = o>,  - co,7,  a 2k  x 8 mHz.  Over  shorter  time  scales,  as  shown  in  the 
inset  of  Fig.  2,  much  of  this  crystal  rotation  takes  place  with  sudden  jumps  in  0crv,  ‘slips’, 
whose  time  scale  is  too  fast  to  be  captured  by  the  top-view  diagnostic.  Intermittent 
behavior  is  a common  feature  in  the  plastic  deformation  (creep)  of  many  materials  (see 
[32]  for  references).  Let  A0m.  denote  the  angular  displacement  between  two  successive 
measurements  of  8<r,..  The  statistics  of  A0fn,  consists  of  a normal  distribution  (from 
measurement  error)  centered  about  zero  with  a width  of  ~ 0.00271,  and  infrequent  larger 
slips.  Due  to  the  known  sign  of  the  i.-beam  torque  and  the  n ambiguity  mentioned  above, 
we  choose  A0C,V  to  lie  in  the  range  [0,7t).  To  separate  statistically  significant  slips  from 
measurement  error  we  further  require  0.00771  < A0rn-  < 0.97ti.  We  find  that  statistically 
significant  slips  account  for  greater  than  90  % of  the  measured  change  in  0lT,,. 

The  I-beam  torque  is  applied  to  all  the  9Be+  ions  in  the  radial  interior  of  the  crystal. 
The  rotating  perturbation,  however,  applies  its  torque  on  the  outer  radial  boundary  of  the 
heavy  ions.  We  therefore  believe  the  stress  due  to  the  competition  between  these  torques 
is  greatest  in  the  region  of  the  heavy  ions  and  anticipate  that  the  slips  of  Fig.  2 are  due 
to  ion  motion  between  the  radial  boundary  of  the  9Be+  ions,  and  the  overall  radial 
boundary  of  the  plasma,  /'o.  This  is  supported  by  the  top-view  images  and  by  simulation 
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FIGURE  2.  Crystal  orientation  0m.  in  the  frame  of  the  rotating  perturbation  for  two  data  sets  with 
different  J_-beam  torques.  The  torque  is  greater  in  data  set  2.  The  lines  are  from  a linear  regression  fit. 
The  inset  shows  a magnified  plot  of  the  first  60  s of  data. 


work  discussed  below,  which  show  most  slips  occurring  as  approximate  rigid  rotations 
of  the  9Be+  ions.  Because  the  simulation  slips  occur  at  a radius  slightly  greater  than 
and  rse  varied  from  data  set  to  data  set,  we  characterize  a slip  amplitude  A„up  by  the 
linear  distance  A $CrvrBe- 

Figure  3 shows  the  distribution  f{Asiip)  of  slips  for  the  two  data  sets  shown  in  Fig. 
2.  Because  we  can  not  distinguish  between  slips  with  amplitude  Asiip  or  Asup  + nnrge, 
where  n is  an  integer,  we  fitted  to  the  function  fs/ip  °c  i^siip  + nKrse)~^  to  determine 
the  agreement  of  the  data  with  a power  law  distribution.  Here  ncu,  is  a cutoff  that  could 
depend  on  the  system  size,  creep  rate,  or  other  factors.  We  obtain  a good  fit  for  any  nCi,t 
but  find  that  y2  g°es  through  a weak  minimum  at  6,  and  therefore  use  that  value 
in  the  following  analysis.  This  results  in  measured  ys  slightly  larger  (<  10  % and  within 
the  uncertainty  of  the  fit)  than  those  resulting  from  nctl,=  1.  In  the  insert  of  Fig.  3 we  plot 
the  measured  power-law  exponent  y as  a function  of  the  creep  rate  Ato  • rge  for  1 0 data 
sets  with  the  same  rotating  perturbation  strength  but  different  _L-beam  torques.  We  find 
that  y decreases  as  the  creep  rate,  a measure  of  the  applied  ±-beam  torque,  increases. 
Decreases  in  the  stick-slip  exponent  with  increased  drive  have  been  observed  in  some 
experimental  systems  [8,  33,  34],  but  not  in  all  [32], 

Inspection  of  Fig.  2 shows  that  the  waiting  periods  (the  time  intervals  between  suc- 
cessive slips)  are  typically  many  seconds.  An  analysis  of  the  waiting  periods  shows  a 
peaked  distribution  with  median  waiting  periods  ranging  from  4 s for  the  highest  ±- 
beam  torques  to  12  s for  the  lowest  ±-beam  torques.  In  Fig.  4 we  plot  the  frequency  of 
slips  versus  waiting  time  for  a data  set  with  moderate  torque.  The  frequency  was  nor- 
malized as  before  by  dividing  the  number  of  slips  nsiip  by  the  binwidth  and  the  duration 
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FIGURE  3.  Distribution  f(Asi,p),  where  f(A,nl,)ilAsnl,  is  the  frequency  of  slips  between  AsuP  and 
Asup  + dAsnp,  for  the  two  data  sets  shown  in  Fig.  2.  The  lines  are  fits  to  a power  law  with  a cutoff  as 
described  in  the  text.  The  insert  shows  the  measured  power-law  exponent  y versus  applied  torque,  as 
parameterized  by  the  creep  rate,  for  all  of  the  data  sets. 


of  the  data  set  (4882  s). 

Earthquakes  are  another  phenomena  that  exhibit  power-law  frequency  versus  ampli- 
tude behavior  and  peaked  waiting  time  distributions.  Although  there  has  been  no  consen- 
sus on  which  is  the  ‘correct’  distribution,  it  has  been  suggested  that  the  waiting-time  dis- 
tributions of  earthquakes  can  be  described  by  lognormal  [35],  Weibull  [36]  or  Gaussian 
[37]  statistics.  In  Fig.  4 we  also  plot  3-parameter  fits  of  the  data  to  these  distributions. 
The  Weibull  distribution,  often  used  to  model  the  time  until  failure  of  brittle  materials, 
is 

ml'"- 1 

p(t)=Ao ^ m > 0,  (2) 

'o 

and  has  a median  waiting  time  value  of  tme,t  = (log2)l'"'/o.  A Weibull  model  for  failure 
can  be  derived  theoretically  as  a form  of  an  extreme-value  distribution,  governing  the 
time  to  failure  of  the  ‘weakest  link’  of  many  competing  processes  [38].  When  the 
shape  parameter  m — 1,  the  Weibull  distribution  reduces  to  the  exponential  distribution 
p(/)  = Aoe~^'^°K  The  exponential  distribution  is  produced  by  a system  with  a constant 
failure  rate  and  thus  no  dependence  on  the  state  of  the  system  (i.e.,  it  is  ‘memoryless’). 
Distributions  where  m 1 result  when  the  outcome  of  a trial  depends  in  part  on  how 
long  the  trial  has  been  conducted. 

As  seen  in  Fig.  4,  the  fit  to  a Weibull  distribution  results  in  a significantly  smaller 
y}  (chi  squared  per  degree  of  freedom),  and  this  is  true  for  the  other  data  sets  as  well. 
If  we  sum  up  the  y}  of  fits  for  all  of  the  10  data  sets,  the  Gaussian  and  lognormal 
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fits  have  statistically  the  same  value  (18.5  and  21.1),  but  the  Weibull  fit  produces 
Xx2  = 10-8.  The  average  and  standard  deviation  of  the  shape  parameter  m is  1.96  ±0.35, 
perhaps  indicating  that  the  ion  crystal  has  a ‘memory.’  For  large  m,  p(/)  approaches  a 
delta  function.  Ward  and  Goes  [36]  have  used  v = 1/m  to  parameterize  the  degree  of 
aperiodicity,  and  found  in  their  earthquake  simulation  data  that  as  they  increased  the 
magnitudes  of  the  quakes  considered,  both  tnied  and  m increases.  We  also  find  these 
trends  in  our  ion-crystal  data. 

Previous  workers  in  the  analysis  of  SGR  ‘starquake’  statistics  have  looked  for  correla- 
tions between  waiting  times  and  the  amplitude  of  the  y-ray  burst.  In  order  to  investigate 
this  with  our  data,  we  sorted  a total  of  918  slips  from  4 data  sets  with  similar  torque  by 
increasing  slip  amplitude,  and  then  into  groups  of  ~ 1 50  slips.  We  next  fitted  the  time 
intervals  of  each  of  these  groups  to  a Weibull  distribution.  We  did  this  both  for  the  time 
interval  until  the  next  slip  of  any  amplitude  (A?+)  and  for  the  interval  from  the  previous 
slip  (At~).  The  results  are  plotted  in  Fig.  5,  and  show  evidence  for  increasing  tmed  for 
A r,  and  much  less  or  no  correlation  for  A t+ . This  is  at  variance  with  results  from  SGRs, 
where  either  no  correlation  has  been  seen  [2 1 , 39, 40],  or  an  anticorrelation  seen  for  A/+ 

[41]- 

Most  experiments  exhibiting  stick-slip  behavior  are  performed  with  ‘constant- 
velocity  driving’  where  the  force  is  applied  through  an  effective  elastic  coupling 
[8,  9,  10].  The  driving  force  of  the  system  is  something  like  F(t)  = K(Vt  — x(t))  where 
‘x’  is  the  ‘position’  of  an  element  in  the  system  (for  example,  the  position  of  a bead  or 
slider  block  in  a chain),  K is  the  effective  spring  constant  that  couples  the  applied  force 
to  each  element  in  the  system,  and  V is  the  constant  average  velocity  that  is  imposed 
on  the  system.  Stick-slip  motion  occurs  for  small  V and  K,  and  a critical  point  exists  in 
the  limit  V— >0  and  K— >0  [15,  33].  If  the  system  gets  stuck,  it  will  eventually  slip  again 
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Slip  Amplitude  (center  of  bin  in  gm) 

FIGURE  5.  Measured  iIIICli  versus  slip  amplitude.  The  lines  indicate  the  mean  value  and  standard 
deviation  measured  for  random  slip  amplitudes. 


because  the  driving  force  increases  linearly  until  slip  occurs. 

However,  our  experiment  is  performed  under  conditions  more  similar  to  ‘constant- 
force  driving’  since  the  J_-beam  radiation-pressure  force  is  constant  in  time  and  applied 
directly  to  the  9Be+  ions.  In  constant-force  driving  the  system  undergoes  a depinning 
transition  at  a critical  force  Fc  and  moves  with  constant  average  velocity  proportional  to 
(F  -Ff)P  for  F > Fc  and  critical  exponent  p [33].  If  a constant-force  system  gets  stuck 
for  a long  period  of  time,  then  this  indicates  F < Fc  and  the  system  should  permanently 
stick. 

In  fact  the  observed  waiting  times  (Fig.  4)  are  long  compared  to  any  known  dynamical 
time  scales  due  to  internal  modes  of  the  system.  Why  doesn’t  the  system  permanently 
stick?  One  possibility  is  that  the  slips  could  be  excited  by  an  external  perturbation. 
By  deliberately  modulating  the  amplitude  of  the  cooling  and  torquing  lasers  we  have 
established  that  the  amplitude  noise  present  in  these  beams  is  not  high  enough  to  trigger 
slips.  Another  source  of  perturbations  is  the  neutral  atoms  of  the  ~ 10-8  Pa  vacuum 
background.  We  have  measured  that  the  relative  abundance  of  9Be+  ions  to  heavy 
impurity  ions  decreases  with  an  e-folding  time  of  ~ 10  hours,  due  to  collisions  with 
background  neutrals.  With  a cloud  of  20,000  ions,  this  gives  a median  time  between 
collisions  of  ~ 1.25  s.  However,  the  relationship  between  these  collisions  and  slips,  if 
any,  is  unclear.  For  example,  if  all  slips  were  caused  simply  by  collisions,  the  waiting 
time  distribution  for  slips  would  be  exponential,  in  contrast  to  what  is  observed. 

We  have  performed  molecular-dynamics  simulations  with  1000  ions  (40  % 9Be+, 
60  % heavy  ions)  with  the  goal  of  better  understanding  the  dynamics  and  origin  of 
the  slips.  About  four  months  of  computation  time  has  been  required  to  simulate  the 
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equivalent  of  ~ 1 .2  s of  experimental  time.  The  simulations  have  produced  three  events 
which  can  be  interpreted  as  slips.  In  these  events,  a rearrangement  of  a small  number  of 
heavy  ions  in  the  vicinity  of  a lattice  defect  produces  a sudden  change  in  the  orientation 
of  the  crystal.  Once  started,  the  slip  eventually  stops  because  the  driving  force  of  the 
_L-beam  is  not  sufficient  to  sustain  continuous  motion.  Because  it  appears  that  the 
rearrangement  of  a few  ions  can  trigger  a slip,  thermal  fluctuations  of  the  ions  may 
be  responsible  for  starting  the  slips  in  both  experiment  and  simulation.  However,  the 
simulations  also  show  the  presence  of  a torsion  mode  of  large  amplitude  and  low 
frequency  (~  10  Hz)  that  produces  a time-dependent  stress  which  may  also  play  a role  in 
triggering  slips.  This  mode  was  not  observed  in  the  experimental  data,  possibly  because 
of  damping  from  laser  cooling. 

In  summary,  we  have  observed  creep  and  stick-slip  motion  in  the  rotational  control 
of  laser-cooled  ion  crystals  in  a Penning  trap.  We  believe  this  system  is  constant-force 
driven  and  may  be  an  experimental  example  of  a subcritical  state  [15]  where  the  slips 
are  triggered  by  thermal  fluctuations  or  by  other  as  yet  unidentified  perturbations  such  as 
collisions  with  neutral  background  atoms.  The  trapped-ion  crystal  system  discussed  here 
possesses  most  of  the  features  of  a self-organized  critical  state  [14, 15].  Therefore  further 
investigations  of  the  stick-slip  behavior  over  a wider  range  of  control  parameters  (_L- 
beam  torque,  temperature,  rotating  perturbation  strength,  background  pressure)  could  be 
useful  for  understanding  the  applicability  of  the  this  concept  to  real  physical  systems. 
Finally,  minimizing  the  occurrence  of  the  slips  is  important  for  some  applications  [5,  6], 
This  can  be  done  by  minimizing  the  J_-beam  torque,  either  through  active  control  of 
the  1-beam  position  or  by  appropriate  tailoring  of  the  1-beam  profile  [42],  Increasing 
the  strength  of  the  rotating  electric  field  perturbation  should  also  decrease  the  frequency 
of  slips  due  to  small  ion  rearrangements.  Two  data  sets  taken  with  half  the  rotating 
perturbation  strength  of  the  data  sets  analyzed  here  showed  an  increase  in  the  number  of 
slips  and  rotational  creep  of  the  ion  crystal. 
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